Low-temperature anomaly in heat capacity 
due to overlapping the spectrums in molecular crystals 
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Investigations of dynamic and thermodynamic properties for a molecular crystal tris-hexafluoro- 
acetylacetonate-iron Fe(02C$IIF§)?, are presented. Heat capacity C P (T) has been measured by 
adiabatic calorimetry method in the temperature range 4.8 — 321 K. An anomaly with a maximum 
at T c = 44.6 K has been discovered. Intermolecular vibrations spectrum was calculated by lattice 
dynamics method in quasiharmonic approximation. Intramolecular frequencies are found by solving 
the Schrodinger equation in approach of small harmonic oscillations. In the frequency interval 
« 30 — 70 cm -1 overlapping the spectrums intra- and intermolecular oscillations has been found. 
The good agreement for calculated and experimental C P (T) occurred to be possible for two sets of 
the force constants. These sets describe two phases above and below T c . The difference between 
phases is connected with freezing of rotation CF3 groups. It has been concluded that the interaction 
between different modes leads to phase transition and anomaly in heat capacity. 
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INTRODUCTION 

Modern materials technology even more often inter- 
feres with study of complicated molecular crystals con- 
taining hundreds atoms in a unit cell. A description of 
molecular crystals is usually based on an idea of inter- 
molecular and intramolecular degrees of freedom. As a 
rule, in simple molecular crystals the spectrums of inter- 
and intramolecular vibrations are energy separated. In 
complicated molecular crystals these spectrums can over- 
lap. Then a possible interactions between different types 
of vibrations cause an occurring of new phenomena (res- 
onances, phase transitions and others). 

In this work the results of investigation of dynamic 
and thermodynamic properties for the molecular crystal 
tris-hexafluoroacetylacetonate-iron Fe(02C§H Fg)^ (or 
Fe(hfac)z, /3-diketonate Fe) are presented. 



EXPERIMENT 

Heat capacity C p (T) of a sample Fe(hf 0,6)3 has been 
measured by adiabatic procedure in the temperature 
range 4.8 — 321 K. In the stepwise heating mode, 141 
experimental points of the heat capacity have been ob- 
tained 0| . After examining the C P (T) dependence, the 
anomalous behavior of heat capacity has been discovered 
in the temperature range 30 — 60 K. It can be clearly 
seen on the plot of derivative dC p (T)/dT, Fig. □ The 
amplitude of an extracted anomaly make up ~ 3% from 
regular heat capacity at T c — 44.6 K (experimental error 
is equal 0.16% within this temperature range). 

To understand the nature of anomaly the static mag- 
netic susceptibility was measured by MPMS-5s SQUID- 
magnetometer (Quantum Design) within the tempera- 
ture range 2 — 300 K. The sample Fe{hf 0,6)3 has 
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FIG. 1: Experimental heat capacity of Fe(hfac)a in coordi- 
nates TdC(T) / dT vs T in the region of anomaly. 



shown paramagnetic behavior, see Fig. [5] Experi- 
mental points coincide well with the Curie- Weiss law, 
x(T) = Ck/(T — Qk), where Curie constant Ck — 
(4.331 ± 0.003) cm 3 Kmol -1 , and paramagnetic temper- 
ature Q K = (-0.1951 ± 0.0002) K. No pronounced mag- 
netic anomalies were observed. One can say, the anomaly 
in heat capacity has nonmagnetic nature. 



INTERMOLECULAR VIBRATIONS 

Computer calculation of thermodynamic properties, 
connected with intermolecular oscillations of Fe(hf 0,6)3, 
was carried out within the framework of lattice dynamics 
method. Earlier this approach has been applied for de- 
scription of dynamic and thermodynamic properties for 
/3-diketonates of metals in the harmonic approximation 
2\. In this paper more complicated procedure based on 
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FIG. 2: Inverse static magnetic susceptibility 1/x as a func- 
tion of temperature. 
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FIG. 3: Intermolecular vibrational density of states at 320 K 
and 10 K. The normalization of curves is made on 6. 



quasiharmonic approach Q has been used. A component 
part of this method is prestress optimization of struc- 
ture. This is caused by lack of information about crystal 
structure and, in the greater degree, about interatomic 
potentials and their parameters. Detailed description of 
optimization procedure and extended possibilities of lat- 
tice dynamics method for molecular crystals had been 
published in Ref. 0. 



The general plan of calculation looks as follows. Initial 
data are the experimental parameters of atom-atomic po- 
tentials, a unit cell and positions of molecules. At first 
the geometry optimization of crystal is carried out. As 
a result, all atoms occupy equilibrium positions at pre- 
determined pressure. Then the dynamic matrix is cal- 
culated. Eigenvalues of this matrix determine oscillation 
frequencies of a crystal and allow one to calculate ther- 
modynamic properties. 

Free energy of the crystal F q h in that approach is calcu- 
lated as F q h = U + F V ib, where U is the potential energy 
and F v n, is vibrational part of free energy 



F mb = - ^2 tiwj (q) + kT ^ In (l 



-hu)j{q)jkT 



hi 



Here k is the Boltzmann constant, uJj(q) are the lattice 
vibrational frequencies. In quasiharmonic approximation 
the free energy for a given crystal structure is taken as 
for harmonic approximation. The anharmonisity causes 
mode frequency dependence of the structural parame- 
ters, making F v u, a function of these parameters as well 
as of temperature. This dependence can be deduced by 
self-consistent manner as a result of optimizing a crystal 
structure while calculating a free energy. 

In this work molecules of a crystal Fe{hfac)^ are con- 
sidered as rigid. The unit cell given in Ref. [5| was used as 
initial structure with parameters a = (9.057 ± 0.004)A, 
b = (13.424 ± 0.005)A, c = (21.591 ± 0.016)A, a = 7 = 
90°, = (116.71 ± 0.02)° for T = 320 K, A unit cell 
contains four molecules Fe{hf 0,0)3. 

For describing the interactions between atoms of dif- 
ferent molecules the Buckingham potential as earlier in 
Ref. was used. In Ref. |(| the universal representa- 
tion of interatomic interaction potentials, Universal Force 
Fields, has been suggested. In terms of Universal Force 
Fields description the Buckingham potential looks as 
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where x — r / p 1 r is the interatomic distance, p is the in- 
teratomic distance at the minimum of the potential (i.e., 
where U{1) = —D), D is the potential well depth and 
£ is a steepness factor. As is conventionally done, the 
general p mn , D mn and Cmn for different m and n types 
of atoms are obtained from homonuclear parameters 
through the use of combination rules p mn = (PmPn) 1 ^ 2 , 
D mn = (D m D n y/ 2 and Cmn = (CmCn) 1/2 - The potential 
parameters used in present work were taken from Ref. [7j 
for light atoms and from Ref. for iron. In the men- 
tioned above representation they are given in the TableHl 
The shape of the cell and the positions of molecules 
were optimized for normal pressure, P — 1 bar, at given 
temperature according to the technique described above. 
These parameters were used as initial ones to calculate 
structural parameters at lower temperatures. Further on, 
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TABLE I: Universal Force Fields parameters of Buckingham 
potential for Fe(hfac)s. 



Atom type, n 


Pn (A) 


D n (kj mol 1 ) 




H 


2.794 


0.127 


11.987 


C 


3.804 


0.332 


14.001 





3.253 


0.512 


13.598 


F 


3.237 


0.218 


14.922 


Fe 


3.957 


1.140 


13.611 




Temperature (K) 

FIG. 4: Heat capacity connected with intramolecular vi- 
brations. Heat capacity without intermolecular contribution 
Cin(T) is shown by circles. Straight-line is heat capacity cal- 
culated on bases of intramolecular vibrations for two Sets of 
force constants. 

the procedure was repeated with stepwise temperature 
decreasing till T = 1 K, Frequencies ujj(q) and density 
of phonon stats (DOS) have been calculated at every tem- 
perature using 5x5x5 points of the Brillouin zone. 

As a result of these calculations, intermolecular vibra- 
tional spectrums and relevant thermodynamic functions 
at different T within the range 1 — 320 K have been re- 
ceived. Intermolecular DOS for 10 K and 320 K is shown 
in Fig. |2| The normalization of curves is made on 6, quan- 
tity intermolecular modes of one molecule. It is seen that 
as the temperature decreases from 320 to 10 K the edge 
of intermolecular vibration spectrum is uniformly shifted 
from » 60 to » 70 cm -1 . Two localized bands separated 



by 6 — 8 cm -1 can be observed near the edge of spec- 
trum. As the temperature varies the marked bands are 
moved. They can consequently take on any values in the 
mentioned frequency range. 

On the base of intermolecular vibrational spectrum 
the corresponding component of heat capacity C mo i (T) 
has been calculated within the range 1 — 320 K . These 
data allow to calculate C,* n (T), the contribution to a 
heat capacity depending from the intramolecular oscil- 
lations. For energy separated intra- and intermolecu- 
lar vibrational spectrums there is the relation Cj n (T) = 
C P {T) — C mo i(T). Calculated in such a way C% n (T) is 
shown in Fig. 0] by circles. 

INTRAMOLECULAR VIBRATIONS 

The vibrational spectrum of a molecule Fe{hfac) 3 is 
found in framework of classical description of small har- 
monic oscillations of atoms near their equilibrium po- 
sitions in molecule. The complex of programs based 
on methods and algorithms published in Ref. [|J has 
been used. The Schrodinger equation for a motion of 
43 atoms with a potential energy U po t = j Uijlilj 
has been solved. Here «,j are the force constants of in- 
teratomic interaction connected with vibrational coordi- 
nates qi and qj. Initial parameters have been structure, 
characteristics of atom positions in the molecule |5J| and 
force constants. Well-known force constants [ToL lllj that 
provide a good agreement of the calculated vibrational 
spectrum with experimental ones were used for valence 
bonds and angles. They were not varied further on. The 
constants connected with non-planar coordinates and 
unknown for Fe{hfac) 3 were varied. The spectroscopic 
data on structural-isomeric compounds of this class as 
well as calculated data for /3 - diketonate chromium 
Cr(02C5H 7 ) 3 were also used. 

In standard method of choosing force constants (based 
on agreement of experimental and calculated spectro- 
scopic data) complications arise when finding force con- 
stants responsible for the low-frequency spectral region 
< 150 — 200 cm -1 . These problems are due to lack of 
any information about frequency spectrum in this region. 
Instead of it, however, it is possible to use experimental 
data on heat capacity in the low-temperature interval 
T < 320 K. In the present work was used new method 
for choosing force constants. The criterion for optimal- 
ity has been the best agreement of functional behavior of 
calculated on the bases of intramolecular vibrations heat 
capacity with Cj n (T). 

As a result of calculation, 123 frequencies responsi- 
ble for interior oscillations of a molecule are found. It 
was discovered, there are two sets (Set 1 and Set 2) of 
force constants gives a good agreement calculated heat 
capacity and Ci n {T). First Set gives a good agreement 
of heat capacities in the temperature range T > T c (see 
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FIG. 5: Low-frequency, q < 210 cm~ , total density of states for two Sets of force constants. Set 1 describes a phase for 
T > T c , Set 2 corresponds to T < T c phase. Intermolecular DOS has 6 frequency modes and calculated for 45 K. Arrows show 
redistribution of intramolecular frequencies for Set 1 and Set 2. 



Fig. f|] Set 1). In the temperature range T < T c a good 
agreement takes place with the second Set (Fig. 0J Set 
2). A disagreement is comparable with experimental er- 
ror. The systematic error above 250 K may comes from 
anharmonic vibrations of intermolecular spectrum. 



DISCUSSION 

Two finding various sets of force constants are ap- 
peared in different functional behavior of a heat capacity 
in areas above and below T c . It gives the basis to suppose 
that these solutions may describe two different structural 
phases. More detailed difference between them can be 
seen on Fig. O Here is shown the total DOS for two Sets 
of force constants in low-frequency region, q < 210 cm -1 . 
One can see that the intramolecular spectrum is sepa- 
rated on parts containing 6 frequencies. It is marked by 
vertical dashed lines in Fig. For Set 1, an overlap- 
ping intermolecular spectrum (6 modes) and intramolec- 
ular spectrum (12 modes) is realized. In T < T c phase, 
Set 2, intermolecular spectrum overlapped only with 6 
intramolecular modes. The principle difference in spec- 
trums takes place only in the region below ss 170 cmT x . 
Just these modes effect essentially on the heat capac- 
ity and low temperature thermodynamic characteristics. 
The vibration modes over rj 600 cm -1 contribute less 
then 1%. The contribution of vibration modes over 



rs 1200 cm is comparable with experimental error of 
CAT). 

More exactly we can say that in phase T <T C (Set 2) 
in spectrums overlapping takes a part an intramolecular 
modes responding for low-frequency bending oscillations 
of (C^Cs-ffi^-plane. Another 6 modes was displaced 
in more high-frequency area and do not take a part in 
overlapping of spectrums. These modes are formed by 
oscillations of CF3 groups. These kind of oscillations may 
proceed to rotational type under increasing of amplitude 
of oscillations. In other words, it is possible that at T < 
T c groups CF3 finish to make rotational movements and 
there is their ordering. 



CONCLUSION 

A new method based on using experimental heat ca- 
pacity data for deciding of force constants allows us to 
take into account a correct density of modes on low fre- 
quency region of spectrum which would be problematic 
on the base of spectroscopic data. We find the over- 
lapping intra- and intermolecular spectrums in complex 
molecular crystal. Spectrums are shifted due to temper- 
ature decreasing and it is possible that due to resonant 
interaction between spectrums arises two different struc- 
tural phases above and below T c . The difference between 
phases is connected with freezing of rotation CF3 groups. 
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And finding out anomaly in heat capacity is the conse- 
quence of spectrums overlapping. 

A theoretical modeling and examining the experimen- 
tal data obtained in this work and earlier |l2( for (3 - 
diketonate chromium lead us to conclusion that overlap- 
ping the spectrums of intra- and intermolecular vibra- 
tions is the characteristic property of molecular crystals 
- diketonate of metals. Certainly our consideration is not 
deprived the lacks occuring from originally separate con- 
sideration intra- and intermolecular fluctuations. More 
perfect technique of simultaneous description of coplex 
crystals is necessary. 
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